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Abstract The aim of this study was to use finite element 
modeling (FEM) as a tool to analyze microwave scattering 
in wood and to verify the model by measurements with a 
microwave scanner. A medical computed tomography scan¬ 
ner was used to measure distribution of density and mois¬ 
ture content in a piece of Scots pine (Pinus sylvestris). 
Dielectric properties were calculated from measured values 
for cross sections from the piece and used in the model. 
Images describing the distribution of the electric field and 
phase shift were obtained from the FEM simulation. The 
model was verified by measurements with a scanner based 
on a microwave sensor. The results show that simulated 
values correspond well to measured values. Furthermore, 
discontinuities in the material caused scattering in both the 
measured and the simulated values. The greater the discon¬ 
tinuity in the material, the greater was the need for compu¬ 
tational power in the simulation. 

Key words Microwave scanning • Wood • Dielectric pro¬ 
perties • Modulated scattering technique • Finite element 
modeling 


Introduction 

Microwave scanning is a fast and nondestructive method of 
measuring internal properties, such as density and moisture 
content (MC) of wood. This method has previously been 
studied, for example, by Johansson et al. 1 The aim of the 
present study was to use finite element modeling (FEM) as 
a tool to analyze microwave scattering in wood and to verify 
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the model by measurements with a microwave scanner. 
Values of phase shift and attenuation from the measure¬ 
ments were compared with the FEM model. 

By transforming the Maxwell equations into second- 
order partial differential equations (PDEs) it is possible to 
solve electromagnetic wave propagation problems. Because 
wood is an inhomogeneous and anisotropic dielectric mate¬ 
rial, the wave scatters in different directions as it propa¬ 
gates. Moisture and heat flow during microwave drying of 
wood have been modeled by Antti, 2 Perre and Turner, 3 and 
Zhao and Turner. 4 No previous models have been made of 
electromagnetic wave propagation inside wood because of 
difficulties in ascertaining the internal structure and be¬ 
cause of the need of computational power. In the present 
analysis, the internal structure of density and moisture con¬ 
tent in wood was determined by computed tomography 
(CT) scanning, as described by Lindgren. 5 From these den¬ 
sity and moisture content values, the distribution of dielec¬ 
tric properties was determined. The finite element modeling 
electromagnetic module in FEMLAB 6 version 3.1 software 
was used to solve PDEs that describe the wave propagation. 


Theory 

The dielectric properties of wood depend on moisture con¬ 
tent, density, frequency, grain angle, and temperature. The 
real and imaginary part of the dielectric permittivity, e' and 
£", perpendicular to the grain versus moisture content and 
density are shown in Figs. 1 and 2. The free water in the cell 
cavities, when above the fiber saturation point (FSP), seems 
to affect permittivity values somewhat differently from 
water that is bound in the cell walls. Parallel to the grain, the 
values are 1.1-2 times higher. 7 

Simulation models were constructed with the finite ele¬ 
ment method as described by, among others, Jin 8 using 
FEMLAB 6 3.1. The models describe how a transverse 
electric (TE) wave propagates through a piece of wood 
surrounded by air. The variations in the material require a 
three-dimensional model to obtain a correct solution. How- 
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Fig. 3. The piece of Scots pine that was used in the study 


Fig. 1. The relative permittivity e' at 10 GHz and room temperature, 
perpendicular to the grain versus moisture content and density 7 



the model and values for e' and e" are interpolated for each 
grid point and time step by the solver. The model uses a 
low-reflectance boundary condition on the transmitting and 
receiving antennas, with electric field E = 1 on the transmit¬ 
ting antenna and with a source field of zero on the receiving 
antenna. The low-reflectance boundary to free space is 
defined as: 

e z -(nxH)+E z =2-E 0z , (3) 

where H is the magnetic field and E 0z = E is the electric 
field. 

Because the incoming wave to the domain is propagated 
parallel to the vertical boundaries and the magnetic field is 
perpendicular to these boundaries, a perfect magnetic con¬ 
ductor (PMC) boundary condition is used. PMC is defined 
as: 


Fig. 2. The relative loss factor e" at 10 GHz and room temperature, 
perpendicular to the grain versus moisture content and density 7 


(nxH) = 0, 

where n is a unit vector normal to the boundary. 



ever, insufficient computational power made it necessary to 
restrict the simulation to a two-dimensional model. Varia¬ 
tions of dielectric properties in the z direction were ne¬ 
glected, which may increase the errors in the final solution. 
In the model domain, the z component of the E field is 
solved by the following equation: 
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where /i r is the relative permeability, which is equal to the 
permeability /r 0 of the free space in the total model domain, 
because wood is not a magnetic material. Furthermore, co is 
the angular frequency, a is the conductivity, and k 0 is the 
wave number in free space. In this case a is close to zero, e* 
is the dielectric permittivity defined as: 

£* = (£'-/•£"), ( 2 ) 

where e' is the relative permittivity and e" is the relative loss 
factor. Tables from Torgovnikov 7 that describe £* at varying 
moisture content, density, and temperature are included in 


Material and methods 

The material chosen for the study was one piece of Scots 
pine ( Pinus sylvestris ) shown in Fig. 3 with dimensions 40 x 
150 x 320 mm. Close to the center of the piece was a cone- 
shaped knot with a diameter varying from 18 to 25 mm. 

A CT scanner (Siemens Somatom) was used to measure 
the density in the cross sections. The geometric shape, di¬ 
electric properties, and moisture distribution in the green 
condition for each cross section were calculated from the 
measured density. 

Density data from CT scanning was used to calculate 
variations in the dielectric properties of the wood and of the 
moisture distribution in the green condition. The difference 
between CT images in green and dry conditions was calcu¬ 
lated and the resulting image was assumed to describe the 
moisture distribution in the cross section. Cross sections of 
the piece, as shown in Fig. 4, were scanned with microwaves 
and CT both parallel and perpendicular to the grain in 
green and dry conditions. 
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Fig. 4. Two cross sections from the object were scanned. One section 
was parallel to the grain and the other section was perpendicular to the 
grain. Regions with low density are darker in the computed tomogra¬ 
phy (CT) images 


Fig. 6. Sensor system for the microwave scanner. The microwave 
signal is modulated by a low-frequency signal at 128 probes in the 
retina. The electromagnetic held at each probe is then extracted from 
the signal. RF, radio frequency; PC , personal computer 
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among others. The measured values of attenuation and 
phase shift in the microwave signal after passing through 
the wood were compared with finite element simulations 
performed in two dimensions. 


Results 


Fig. 5. Brief description of the working procedure in the generation of 
finite element models 


Before the moisture content was calculated from the 
CT images, it was necessary to compensate the image of 
dry wood for shrinkage and deformation. Transformation 
of the CT image of dry wood to the shape of green wood 
was done by using elastic registration as described by 
Sorzano et al. 9 Figure 5 shows how the FEM model was 
generated using density images in green and dry conditions 
when the electric field (E-field) was oriented perpendicular 
to the grain (Fig. 4). The same procedure was used when the 
E-field was parallel to the grain, but the values for e* were 
multiplied with a moisture-dependent factor. 7 Within the 
knot, the E-field was perpendicular to the grain in both 
models. 

A microwave scanning system described by Johansson 10 
that measures attenuation and phase shift of the transmitted 
electromagnetic field every 8 mm was used for the measure¬ 
ments. Figure 6 shows a schematic drawing on the main part 
of the system. The wood is illuminated by a quasiplane wave 
generated by a slotted waveguide that acts as a transmitting 
antenna where the electrical field is perpendicular to the 
sensor (9.375 GHz). 11 The wave is locally perturbed at di¬ 
poles in the retina by a low frequency signal. This will cause 
a modulation of the radio frequency (RF) signal that is 
proportional to the electromagnetic field at the dipoles. 
This technique, known as modulated scattering technique 
(MST), has been described by Bolomey and Gardiol 12 


The results show that simulated values correspond well to 
measured values. Figures 7-10 show attenuation and phase 
angle in the simulated and measured field after transmit¬ 
tance through dry wood. The measured values were ob¬ 
tained by repeated measurements with increasing antenna 
distance. Figures 11 and 12 show the attenuation when the 
E-field is oriented perpendicular and parallel to the grain. A 
clear pattern from the knot in the phase shift for dry wood 
is shown in Figs. 13 and 14. The phase shift perpendicular to 
the grain for green wood in Fig. 15 also shows a pattern 
caused by the knot. The knot is not visible in simulated 
or measured data for green wood parallel to the grain 
(Fig. 16). 


Discussion 

The phase shift is periodic and can only be measured and 
simulated in the interval (-71, n). The vectors were un¬ 
wrapped by changing absolute jumps greater than ;rto their 
2 n complement. This only works when the phase shift be¬ 
tween adjacent points is less than 2k. The measurements 
that deviated from the model were attenuation in green 
wood with the E-field perpendicular to the grain and phase 
shift in dry wood with the E-field parallel to the grain. The 
first case can be explained by variations in the z direction. 
The moisture content was higher on both sides of the 
scanned cross section. In the second case, the result is very 
sensitive to changes in the correction factor for permittivity 
when the E-field is parallel to the grain. 






































































409 


\ Iransm tiling antcima I 



rcccr.mg anieima 


Fig. 7. Simulated attenuation of the electric field (E-field) in dB 
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Fig. 8. Measured attenuation of the transmitted E-field 



Fig. 9. Simulated phase angle of the E-field in radians 


Fig. 11. Simulated and measured attenuation of electromagnetic (EM) 
wave with the E-field oriented perpendicular to the grain after trans¬ 
mittance through green and dry wood 
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Fig. 12. Simulated and measured attenuation of EM wave with the E- 
field oriented parallel to the grain after transmittance through green 
and dry wood 
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Fig. 10. Measured phase angle of the transmitted E-field 



Conclusions 

The model corresponds well to the measured values, which 
show that the FEM simulation could be a useful tool for 
analyzing microwave scattering in wood. Because wood is 
an inhomogeneous material, there are discontinuities in 
density and moisture content in the wood pieces that result 


Fig. 13. Simulated and measured phase shift of EM wave with the E- 
field oriented perpendicular to the grain after transmittance through 
dry wood 

in a huge scattering in the measured and simulated values. 
The greater the discontinuities, the greater is the need for 
computational power in the simulation. Variations in the z 
direction give errors if the models are restricted to the xy 
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Fig. 14. Simulated and measured phase shift of EM wave with the E- 
held oriented parallel to the grain after transmittance through dry 
wood 



Fig. 15. Simulated and measured phase shift of EM wave with the E- 
held oriented perpendicular to the grain after transmittance through 
green wood 


plane. The model for dry wood is sensitive to changes in the 
correction factor for permittivity. 



Fig. 16. Simulated and measured phase shift of EM wave with the E- 
held oriented parallel to the grain after transmittance through green 
wood 
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